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Time-reversal invariant three-dimensional topological insulators can be defined fundamentally by 
a topological field theory with a quantized axion angle of or tt. It was recently shown that 
fractional quantized values of 9 are consistent with time-reversal invariance if deconfined, gapped, 
fractionally charged bulk excitations appear in the low-energy spectrum due to strong correlation 
effects, leading to the concept of a fractional topological insulator. These fractionally charged 
excitations are coupled to statistical gauge fields which ensure that the microscopic degrees of 
freedom, the original electrons, are gauge-invariant objects. A first step towards the construction of 
microscopic models of fractional topological insulators is to understand the nature of these statistical 
gauge theories and their corresponding phases. In this work, we show that low-energy effective gauge 
theories of both Abelian or non-Abelian type are consistent with a fractional quantized axion angle 
if they admit a deconfined phase or a Higgs phase with gauge group broken down to a discrete 
subgroup. The deconfined phases support gapless but electrically neutral bulk excitations while 
the Higgs phases are fully gapped. The Higgs and non-Abelian deconfined phases exhibit multiple 
ground states on boundaryless spatial 3-manifolds with nontrivial first homology, while the Abelian 
deconfined phase has a unique ground state. The ground state degeneracy receives an additional 
contribution on manifolds with boundary due to the induced boundary Chern-Simons term. 

PACS numbers: 73.43.-f, 75.80.-|-q, 71.27.-|-a, 11.15.-q 



I. INTRODUCTION 

Topological insulatorsir~ are new states of quantum 
matter that cannot be adiabatically connected to conven- 
tional insulators. They are fully gapped in the bulk but 
support gapless boundary modes which arc protected by 
discrete symmetries Topological insulators were first 
discovered in HgTe quantum wells J^'-iS More recently, 
three-dimensional (3D) time-reversal invariant topolog- 
ical insulators have attracted a great deal of attention, 
for the most part due to the theoretical predictio n^'^d^ 
and subsequent experimental detectiorti^^— of their pro- 
tected helical surface states. The existence and stability 
of these surface states is protected by a bulk Z2 topo- 
logical invariant which corresponds physically to a topo- 
logical magnetoelectric effect^ (TME) . The experimental 
observation of the TME, for instance by way of magne- 
tooptical measurements is a key goal in the field 
which is being actively pursuedi21"i2£ From a theoreti- 
cal standpoint)^ the TME is described at energies much 
smaller than the energy gap by the addition of a term 
proportional to E • B (with E and B the electric and 
magnetic fields, respectively) to the usual Lagrangian for 
Maxwell electrodynamics, i.e. axion electrodynamics*^ 
This term is in fact the Abelian version of the topologi- 
cal 0-term in quantum chromodynamicS )^^'^^ and its co- 
efficient is periodic (under certain conditions^) with 
period 27r. Since E • B is odd under time-reversal sym- 
metry (T), the only values of 6 allowed by T are or 
TT mod 27r, with ^ = for the Z2 trivial insulator and 
6 = TT for the Z2 nontrivial insulatorj^ Topological insula- 
tors can be described microscopically by noninteracting, 



spin-orbit coupled electrons hopping on a lattice, and the 
axion angle 9 can be computed from a knowledge of the 
single-particle wave functions in momentum space^i^ 

In the presence of electron interactions, the concept of 
a topological insulator can no longer be defined in terms 
of band theory, but the fundamental definition in terms 
of the topological field theory and the axion response 
function remains generally validi^ Weak interactions can 
either turn the noninteracting helical surface state into a 
weakly interacting helical Fermi liquid with spin-charge 
coupled collective modes^^ or drive a transition to a 
superconducting statei^i^ Strong enough interactions 
can lead to spontaneous T-breaking on the surface and 
ferromagnetic4i2^"— or helical spin density wave order 4^ 
In contrast, the bulk of a topological insulator is fully 
gapped and thus expected to be perturbatively stable to 
interactions. On the other hand, exotic states known 
as topological Mott insulators^"— have been theoreti- 
cally proposed, whereas a topologically nontrivial band- 
structure is dynamically generated as a consequence of 
strong electron-electron interactions. Although these are 
strongly interacting states, their mean-field description is 
still that of a topological band insulator, and the axion 
angle 9 remains quantizedi^*' to or tt mod 27r. Another 
type of topological Mott insulator has been theoretically 
propose d^^'^^ in which spin-charge separation leads to a 
bulk insulator with a helical liquid of gapless spinous, but 
an electromagnetic 0-term is not generated because the 
spinous are electrically neutral. 

In many regards, the 3D T-invariant topological insu- 
lator can be viewed as a generalization of the 2D integer 
quantum Hall effect (IQHE) to 3D. The topological Z2 
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quantization of the bulk axion angle 9 in 3D is the di- 
rect analog of the topological Z quantization of the bulk 
Hall conductance in 2D. By analogy with the relation be- 
tween the IQHE and the fractional quantum Hall effect 
(FQHE), one is naturally led to the question whether 
there can exist a "fractional 3D topological insulator" 
which preserves T but is characterized by a fractional 
quantized axion angle, i.e. where is a non-integer, ra- 
tional multiple of tt. A T-invariant fractional topologi- 
cal insulator was first proposed in two dimensions, and 
more recent works have shown the robustness of such a 
statei^"— In the special case of conserved Sz, the QSH 
insulator is equivalent to two decoupled IQHE systems 
with equal and opposite effective magnetic fields. Both 
IQHE systems can be driven into FQHE states by adding 
electron-electron interactions and appropriately tuning 
the effective magnetic fields while keeping them equal 
and opposite, which yields a fractional QSH insulator 
without breaking T . On the other hand, this procedure 
is not directly applicable to three dimensions, where one 
cannot in general reduce a T-invariant topological insu- 
lator to two decoupled topological states which break T 
in an equal and opposite wayi^ 

In our previous workj^ we introduced a theory of 
T-invariant fractional topological insulators (fTI) in 3D 
based on a parton constructionj ^^'^^ Simply postulating 
that an electron, under the influence of strong interac- 
tions in the underlying lattice Hamiltonian, fractional- 
izes into "colors" of partons gives a realization of a 
fractional topological insulator, as long as each color of 
partons forms a topological band insulator. In order to 
ensure that outside the fTI the partons recombine into 
electrons, we needed to introduce additional degrees of 
freedom, analogous to the statistical gauge fields in the 
theory of the FQHE.^^ Outside the topological insulator 
these gauge fields are confining. As a consequence, the 
partons, which are charged under the "statistical" gauge 
fields, can never be observed in isolation but are instead 
confined into electrons, just like quarks are confined into 
mesons and baryons in quantum chromodynamics. The 
theory of the partons together with the statistical gauge 
fields should be thought of as a low-energy effective de- 
scription of the system. The goal of this paper is to 
show that consistent low-energy theories with the char- 
acteristic properties of a fTI exist and to analyze their 
properties]^ 

There are many choices and questions associated with 
the additional statistical gauge fields. Are they Abelian 
or non- Abelian in nature? By construction they should 
confine outside the fTI, but what phase do they realize 
inside the fTI? In Ref. [SJ we mostly focused on the case 
of a deconfined gauge field. It is easy to write down de- 
confined models and to demonstrate the theoretical con- 
sistency of fTIs in principle. Such models can either be 
Abelian or non-Abelian. Examples of both types were 
presented in our earlier work. One generic feature that 
all deconfined realizations of fTIs share is the presence of 
additional gapless degrees of freedom. One typically con- 



siders topological insulator phases that are fully gapped 
in the bulk. In the case of a deconfined realization of 
a fTI we have to slightly generalize this understanding 
by demanding that all degrees of freedom charged under 
electromagnetism are gapped, while allowing for electri- 
cally neutral gapless degrees of freedom. In that case, the 
system is indeed an insulator as far as electrical transport 
is concerned. These additional neutral gapless degrees of 
freedom play a role similar to the low-energy phonons in 
a band insulator, and do not spoil the quantization of the 
electromagnetic response. 

An alternative to the deconfined realization of the fTI 
is to put the statistical gauge fields in a Higgs phase. We 
will present an explicit model of this kind in this work. 
In this case, the system is truly gapped with no gapless 
degrees of freedom. The statistical gauge group is broken 
down to a discrete gauge group, and the latter is sufficient 
to ensure that in physical states the net electron number 
is integer. 

As we will elaborate in more detail in this work, the 
basic topological feature of the fTI, a fractional an- 
gle, is completely robust in that it only depends on the 
total number of partons and not on any of the details 
of the statistical gauge sector — whether it is Abelian 
or non-Abelian, deconfined or Higgsed. This universal 
9 angle characterizes the TME in the material together 
with its physically measurable surface properties, such as 
a surface FQHE with half the conductivity of a typical 
Laughlin state per surface as well as the corresponding 
magnetooptical Kerr and Faraday effects. 

There is, however, a second topological feature of the 
fTI that does depend on details of the gauge field sector: 
the ground state degeneracy. Some basic features of the 
ground state degeneracy have already been discussed in 
Ref. 54, but here we will elaborate on this. It has re- 
cently been proven^ that for a fully gapped system, 9 
can only be fractional if the ground state on the 3-torus 
is degenerate. This observation makes it clear that a 
confined phase is not an option for the fTI, as it would 
result in a unique ground state in a completely gapped 
system. In the case of a deconfined realization of the fTI, 
we will find that the ground state in the Abelian mod- 
els is unique, whereas in the non-Abelian models it is 
typically degenerate. The reason why this is not in con- 
tradiction with the theorem of Ref. [5^ can be traced to 
the fact that we have additional gapless degrees of free- 
dom, therefore violating the assumption of Ref. HI that 
the system is gapped. We will demonstrate ex plic itly 
how this allows us to avoid the arguments of Ref. (591. In 
the Higgs models, the ground state is degenerate as the 
aforementioned theorem requires, and we determine the 
ground state degeneracy. 

We illustrate our findings with the study of three ex- 
amples. All three models correspond to A^c = 3 partons. 
Model A is an Abelian, deconfined realization of a fTI 
based on a.U{\) x U{1) statistical gauge field. Model B 
is a non-Abelian, deconfined realization based on S'[/(3). 
Model C is a Higgsed realization where the gauge group 
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of model B is broken down to Z3 . 

The structure of the paper is as follows. In Sec. [TTl we 
introduce three general classes of parton-gauge boson ef- 
fective theories: the Abelian U{1) deconfined models, the 
non-Abelian deconfined models, and the Higgs models. 
We illustrate the general ideas in each case by a specific 
example of three colors of partons {Nc = 3), correspond- 
ing to three distinct classes of fTI but all with a fractional 
axion angle 9 = 7r/3. In Sec. lIIIi we prove the topological 
quantization of the axion angle in fractional multiples 
of 77. We derive the quantization of the axion angle in 
both Abelian and non-Abelian models, whether gapless 
or gapped. In Sec.lIV| we discuss the issue of ground state 
degeneracy on spatial 3-manifolds of nontrivial topology. 
This is the direct analog of the topological degeneracy in 
FQHE states on Riemann surfaces,^" and can be taken as 
an indication that fTI states exhibit topological order in 
the many-body sensci^ In Sec.|V]we briefly speculate on 
the nature of the gapless surface states, and summarize 
the paper in Sec. I VII Appendix A and B clarify some 
technical issues concerning the ground state degeneracy 
on a 3-manifold with boundary. 

II. THREE CLASSES OF PARTON MODELS 

In this section, we introduce the three basic classes of 
possible statistical gauge sectors for an fTI, and describe 
in more detail a specific Nc = 3 model in each class. 

A. Deconfined models 

As described in the Introduction, one way to realize a 
fTI is to drive the statistical gauge fields into a decon- 
fined phase. This can be realized either in an Abelian 
or a non-Abelian setting with rather different properties. 
What is common to any deconfined realization of a fTI is 
the appearance of extra gapless matter. This extra gap- 
less matter is neutral from the point of view of the f/(l)om 
Maxwell gauge field. Therefore, the system is still an in- 
sulator, i.e. all degrees of freedom charged under C7(l)cm 
are gapped. These additional gapless degrees of freedom 
can be considered as soft "phonons", since they do not 
enter the electrically charged sector. 

1. Abelian models 

Model A: In the simplest Abelian model of a fTI with 
Nc = 3, the electron fractionalizes into three fermionic 
partons ipi, i — 1,2,3 with a statistical U (1)^ ~ U {1)a x 
U{1)b gauge group in addition to the C/(l)om Maxwell 
gauge group. The charge assignments are given in Ta- 
ble U where e is the electromagnetic gauge charge and 
g is the f/(l)^ "statistical" gauge charge. These gauge 
groups ensure that the only gauge-invariant operator that 
carries Maxwell charge is the product of the three parton 
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TABLE L Gauge charge assignments of parton fields un- 
der electromagnetic ?7(l)em and statistical U{1)^ = U{1)a x 
U{1)b gauge groups in the simplest Abelian model for Nc — 3. 



operators, hence the gauge-invariant electron operator is 
''/'iV'2'03- Indeed, the generators of the weight lattice (see 
Sec. IIIIB 2p corresponding to the representation of Ta- 
ble U are ei = {2g,-g), 62 = {-g,2g), 63 = {-g,-g), 
and the equation niG; — requiring an operator 
'0i^'02^V'3^ to be gauge- invariant has the one-parameter 
family of solutions ni ^ n2 = n^, i.e. operators of 
the form ('!/'i'02'03)" of which only n — 1 corresponds 
to an operator with electromagnetic charge e. For a T- 
invariant system, all three partons have a real mass. As 
we will review in the next section, in the topologically 
nontrivial phase all three partons have a real but nega- 
tive mass. 

This model is the simplest of the U{1)^'' /U{l)dia.g 
Abelian models described in Ref. [HJ. The generic model 
starts out with a statistical U{1)^'' gauge group and 
Nc partons. Each parton carries charge e/Nc under 
the Maxwell ?7(l)om gauge group. The ith parton car- 
ries charge g under the ith U{1) factor of the statistical 
gauge group, and is neutral under the remaining U{1) fac- 
tors. One then takes the quotient of the statistical gauge 
group by its diagonal t/(l)diag subgroup whose genera- 
tor is simply the sum of the generators of the individual 
U{l)i factors. For the particular case of Nc — 3, the 
gauge group presented above represents the two remain- 
ing gauge group factors U{1)a = 2U (l)i - C/(l)2 - U{1)3 
and;7(l)B = -U{l)i+2U{1)2-U{l)3. These are two lin- 
ear combinations that are orthogonal to t/(l)diagj in the 
sense that their generators satisfy trt^iiom = triB^Gm = 
0, with tA = dia,g{2g, -g,-g), Ib = d\a.g{-g,2g,-g), 
and iom = diag(e/3, e/3, e/3), and tr is the trace in 
the representation of Table [J For this choice of gen- 
erators, one obtains the convenient feature that the sta- 
tistical parton charges are integer multiples of g. These 
generators are however neither orthogonal to each other 
(tri^is ^ 0) nor properly normalized {trtA,BtA,B 7^ 9^)- 
Therefore, if one really starts out with a gauge 
theory, the Maxwell gauge kinetic term trF^^F^^'^ would 
contain a mixed F^^Fg'^ term. The generators can how- 
ever be easily orthonormalized, and upon doing so we find 
precisely the generators Hi and H2 [Eq. (|13p ] of the max- 
imal diagonal subgroup U{1)'^ of S';7(3) (see Sec. lIIIB2p . 
The same construction works for general number of "col- 
ors" Nc, yielding once again parton charges which can 
be chosen to only take values 2g and —g under the vari- 
ous Abelian factors, at the price of non-diagonal kinetic 
terms. As the gauge kinetic terms do not impact the 
topological properties this is of no importance. This al- 
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lows one to study a simplified, topologically equivalent 
version of the model where one takes the above charge as- 
signments with standard Maxwell terms (no mixing) and 
forgets about the fact that the U{1)^''~^ gauge group 
originated from J7(l)^'= in the first place. If desired, a 
set of orthonormal generators is provided by the Cartan 
generators Hi, . . . ,Hm^-i of SU{Nc), i.e. the genera- 
tors of its maximal diagonal subgroup J7(l)^'=~^. In fact, 
as far as the electrically charged degrees of freedom are 
concerned, the J7(l)^'=/J7(l)diag model discussed above 
is equivalent to a Higgs SU{Nc) model where SU{Nc) is 
spontaneously broken to its maximal diagonal subgroup 
f/(l)^--i. 

Upon integrating out the massive partons and setting 
the electromagnetic gauge potential to zero, in the fTI 
phase we are left with a pure U{1)^''^^ dynamical gauge 
theory with 6 = t: mod 27r. For a continuum theory, 
there is a single phase, the deconfincd Coulomb phase. 
This phase corresponds to a free field infrared fixed point 
at which the renormalized coupling g vanishes. In other 
words, in the infrared we have free Dirac fermions (the 
partons) and photons (the statistical gauge fields). For a 
theory defined on the lattice, we obtain two phases, the 
deconfined Coulomb phase and the confined phasej^"— 
The fTI phase corresponds to the Coulomb phase, which 
is equivalent to the continuum theory in the infrared ex- 
cept for a doubling of the number of fermion species. The 
topological properties are the same in both the contin- 
uum and lattice cases. However, the value (0 or tt mod 
2tt) of the 9 angle for the (free) partons is dependent upon 
the choice of regularization procedure in the continuum 
theory, while it is fixed by the lattice Hamiltonian in the 
lattice theory.— In analogy with the results of Cardy and 
Rabinovici^ and Cardy^ for Zjv gauge theory in 3 + 1 
dimensions, one may however wonder whether the statis- 
tical O-teim with nonzero 9 = tt could affect the phase 
diagram of the statistical gauge sector, since the gauge 
theory is now dynamical. This is not so because a U{1) 
gauge theory can be viewed as the N oo limit of a Zjv 
gauge theory, and in this limit the "electric" charges of 
Zjv gauge theory disappear from the spectrum^ and the 
0-term has no effect on the bulk free energy. Therefore, 
the Abelian models are simplest to analyze. For exam- 
ple, the free fields do not allow for a nontrivial ground 
state degeneracy. 



2. Non-Abelian models 

One may question whether a deconfined phase is actu- 
ally realizable for non-Abelian models. Although pure 
SU{Nc) Yang-Mills theory is generally believed to be 
confining at zero temperature, adding gapless matter can 
result in a deconfined phase. For example, Af — 4: su- 
persymmetric Yang-Mills (SYM) theory with any gauge 
group, in particular with gauge group SU{Nc), is ac- 
tually conformally invariant and hence in a deconfined 
phase. In fact, this is the generic case even in non- 



supersymmetric Yang-Mills theories with enough gapless 
matter: for matter in the fundamental representation as 
well as for matter in two- index tensor representations it is 
believed that SU{m), SO{m) and Sp{m) gauge theories 
all exhibit a "conformal window" , that is, the gauge the- 
ory is in a conformally invariant and deconfined phase as 
long as the number of flavors Nf is within a certain finite 
range. A conjecture for the exact values of Nf that bound 
the conformal windows has been put forward for example 
in Ref . 1691. The upper end of the conformal window is the- 
oretically well established: for Nf larger than this maxi- 
mal value the gauge theory loses asymptotic freedom. It 
is weakly coupled in the infrared but ill-defined in the 
ultraviolet. Just below the upper end of the conformal 
window it can be established, using perturbation theory, 
that at least for large Nc the theory indeed flows to a con- 
formal fixed point.— The lower bound of the conformal 
window is mainly conjectural based on partial resumma- 
tions. There has been a lot of recent activity on numer- 
ical studies of the conformal window using lattice gauge 
theory. While the precise lower bound of the conformal 
window is still up to debate, the existence of a conformal 
window has been firmly established. In gauge theories 
with M = 1 supersymmetry the full conformal window 
has been mapped out using the power of holomorphyf^^ 
The J\f — 1 supersymmetric SU{Nc) gauge theory is for 
example conformal for |iVc < Nf < 3Nc. 

As for the Abelian models, the deconfined non-Abelian 
models have additional gapless degrees of freedom. For 
example, in the case of A/" = 4 SYM theory the extra mat- 
ter consists of adjoint fermions and scalars, in addition 
to the gauge fields. One important difference compared 
to the Abelian case is that in the case of a non-Abelian 
gauge theory at a nontrivial fixed point the extra matter 
is not free, but remains interacting with a fixed renor- 
malized coupling , where 5* is fixed by the requirement 
that the renormalization group /3-function vanishes. In 
the special case of A/" = 4 SYM theory, g^, is in fact a free 
parameter and the /3-function vanishes identically for all 
values of the coupling, -- Not only do the statistical gauge 
fields remain strongly coupled to each other, the partons 
remain strongly coupled to the statistical gauge fields. 
Below we will argue that the value of the fractional 9 an- 
gle only depends on the number of partons the electron 
fractionalizes into, and is completely robust even against 
these strong interactions that remain in the low-energy 
effective theory of the partons. In fact, in case of A/" = 4 
SYM theory this has been demonstrated explicitly in the 
extreme limit of very large coupling, employing a holo- 
graphic realization of this particular model,-'^ 

Model B: The canonical example in this class, which 
we will refer to as an example in various sections, is a 
model with Nc = i partons of electric charge e/3 cou- 
pled in the fundamental representation to a SU{3) gauge 
field with some additional massless adjoint matter fields 
to drive it into the deconfincd phase. Details of this ex- 
tra matter are largely irrelevant from a topological point 
of view, except for the question of the ground state de- 
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generacy on T^. M = A SYM theory can serve as the 
canonical example. (It should be clarified that the exis- 
tence of massless matter fields is consistent with the par- 
tons being gapped. The massive partons, which carry the 
fundamental representation of the SU{3) group, are the 
only electrically charged particles and dominate electro- 
magnetic transport. Being gapped, they are irrelevant 
for the low-energy dynamics of the statistical gauge field. 
The latter is driven by the massless SU{3) charged but 
electrically neutral extra matter.) 

One consequence of the residual interactions in the in- 
frared is that sometimes T is spontaneously broken in the 
topologically nontrivial phase, that is at = tt. Whether 
this happens or not depends on many details of the sta- 
tistical sector.''^ For the purpose of constructing parton 
models of fTI, the important lesson to remember is that, 
even though spontaneous T-breaking does occur some- 
times, non-Abelian gauge theories with unbroken T do 
exist. In the large Nc limit the situation is slightly bet- 
ter understoodi^^ while pure SU{Nc) Yang-Mills theory 
is believed to spontaneously break T at = tt, additional 
gapless matter can prevent this. In particular, = 4 
SU{Nc) SYM theory is T- invariant at 6* = tt in the large 
Nc limit. 

One other important difference between the Abelian 
and the non-Abclian models concerns spin. In model 
A, each of the three partons belongs to a different one- 
dimensional irreducible representation of the statistical 
gauge group C/(l)^. The requirement that for a fermionic 
state the many-body wave function has to be antisym- 
metric under interchange of both color and spin indices 
does not constrain the symmetrization properties of the 
spin quantum numbers, and hence the total spin, because 
gauge invariance does not require a separate, complete 
antisymmetrization of the color indices in this case. A 
spin-1/2 electron is always possible. In model B however, 
gauge invariance under the SU{3) group requires that the 
many-body wave function be completely antisymmetric 
in color indices. Correspondingly the spin indices have to 
be completely symmetrized, and in model B the electron 
would have spin 3/2. Because spin rotation invariance is 
already broken in a topological band insulator by spin- 
orbit coupling, this does not create a problem, but it is 
certainly an aspect of our non-Abelian models to keep in 
mind. 



3. General deconfined model 

The general deconfined model will have some Abelian 
(free) factors and some non-Abelian (interacting) factors. 
A large class of models of this type was introduced in our 
earlier workiH The statistical gauge group in this gen- 
eral model is n^ii )/C^(l)diag- The total number 
of partons is 

Nc^Y. ^c- (1) 

a=l 



For every flavor a = l,...,Nf we have A''" partons 
transforming in the fundamental representation of the 
non-Abelian SU{Nc) factor and carrying [/(1)° charge 
g° — g/N^ as well as electromagnetic charge q°™, while 
being neutral under all the other gauge groups. The di- 
agonal subgroup C/(l)diag, whose generator is the sum of 
all the ?7(1)° generators, is modded out. The only gauge- 
invariant operator that carries nonvanishing electromag- 
netic charge is the product of all the partons. To be 
gauge- invariant under each of the individual non-Abelian 
factors, one needs to form baryonic operators out of the 
partons in the fundamental representation of that fac- 
tor. These individual baryons however will carry U{l)a 
charge g. Only the product of all the individual baryons 
is gauge- invariant, as this is the only way to get a gauge- 
invariant operator whose only statistical U{1) charge is 
the charge under the diagonal subgroup J7(l)om, which 
is removed from the statistical gauge group. The electro- 
magnetic charge of this gauge invariant operator is 

Q^E^c^r, (2) 

o=l 

hence the g®™ have to be chosen in such a way that Q = e. 

Another alternative one might want to consider is to 
use models based on orthogonal or symplectic gauge 
groups. The gauge- invariant operators in such theories 
are the mesons QaQbS'^'' and qaqbi^""^, respectively, where 
S"''' is the Kronecker delta and cj"'' the corresponding an- 
tisymmetric invariant tensor for the symplectic group. 
Both theories have baryons Qai ■ ■ ■ 9ajv e"i" "™<: . However, 
in the symplectic theory Sp{Nc) where Nc is even, these 
baryons are not independent operators but are equivalent 
to a product of mesons, because the symplectic tensor w"'' 
is antisymmetric. By analogy with the [U{M) x Sp{2k)]i 
Chern-Simons theory of the parafermion FQH states 
obtained from a parton construction;^ these theories 
hold the promise of generating more exotic surface states. 
However, if we take the partons to have electric charge q, 
not only we obtain baryons of electric charge NcQ which 
would suggest q = e/Nc as before, but the mesons are 
also charged as we can make mesons from two fundamen- 
tal quarks. In the SU{Nc) case, mesons are made from 
a fundamental quark and an antifundamental antiquark. 
Therefore, this time we are forced to identify the mesons 
with the electrons J- Furthermore, we have to assign the 
partons charge g/2, and at the same time take the num- 
ber of partons Nc to be even so that the baryons carry 
an integer multiple of the electron charge, which is con- 
sistent with the fact that symplectic groups are only de- 
fined for even Nc- To ensure that this mesonic electron is 
a fermion, we need the electron to split into two different 
partons, a fermion and a boson, both in the fundamental 
representation of the gauge group. We would also get 
gauge-invariant scalar boson with the same charge e as 
the electron, but presumably these can be made gapped. 
In this case, the axion angle 9 has to be an integer mul- 
tiple of NcTt/A from the Nc fermions of charge e/2. This 
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TABLE II: Gauge charge assignments of parton fields un- 
der electromagnetic (7(l)cm and unbroken statistical Za gauge 
groups in the simplest Higgs model for Nc = 3. The third col- 
umn corresponds to the triality of the representation. 



structure is somewhat reminiscent of the Z2 spin liquid 
model of a fTI put forward in Ref. [sol . 



B. Higgs models 

In order to get a completely gapped system that re- 
alizes a topological insulator, it is sufficient to add elec- 
trically neutral Higgs fields to the deconfined models of 
Sec. Ill Al and consider a specific pattern of spontaneous 
symmetry breaking G — >■ iJ with G the original gauge 
group, such that the unbroken gauge group H in the 
Higgs phase is discrete. One can modify model B above 
to realize this possibility: 

Model C: A simple Higgs model can be obtained 
by augmenting the simplest non-Abelian three-parton 
model, model B above, by inclusion of two complex scalar 
Higgs fields in the adjoint representation of SU{3). 

The center of SU{3) is Z3 and the adjoint Higgs fields 
are neutral under this center symmetry. If one of the 
scalars acquires a vacuum expectation value, the contin- 
uous part of the gauge group is broken to its maximal 
diagonal subgroup and one recovers the fermions 

and gauge fields of model A together with several charged 
scalars. Giving generic noncommuting expectation val- 
ues to both adjoint Higgs scalars completely breaks the 
continuous part of the gauge group and generically only 
the discrete Z3 center symmetry is unbroken. Under 
this discrete gauge group the parton charges are given 
in Table III Bl The third column corresponds to the 
nonzero triality fc = 1 mod 3 of the fundamental rep- 
resentation with character e'^'^^^^. In general, for a 
SU{Nc) — Zjv^ Higgs mechanism we want the partons to 
be in a representation of SU (Nc) with nonzero A'c-ality 
k = l,2,...,Nc-l mod with character 6^'^^'=/^- ^ 1. 
This discrete subgroup of the original continuous statisti- 
cal gauge group is completely sufficient to ensure that all 
gauge-invariant operators have charges that are an inte- 
ger multiple of the electron charge. Indeed, the center Z3 
of SU{3) is generated by a single element of the Cartan 
subalgebra, H2 oc diag(l,f,— 2) in Eq. ([T^. The statis- 
tical weight lattice is generated by the (unnormalized) 
weight vectors ei — 1, 62 — 1, 63 = —2. A statistically 
neutral operator -0"^?/'2 ^-i/^g^ must satisfy X^i'^i^i = 
which leads to ni + n2 = 2n3 , hence the total Maxwell 
charge is e(ni + n2 + '^3)/3 = n^e € Ze. V'i'02V'3 is still 
the simplest gauge-invariant operator, but as long as the 



fermions have internal spin states that can be antisym- 
metrized, ipf or ipf'^'j with i,j = 1, 2, 3 and i ^ j would 
give rise to additional gauge-invariant operators. Phys- 
ically these operators can be considered as an electron 
combined with condensed Higgs bosons. In the relativis- 
tic continuum theory the fermions carry spin 1/2 and 
■01 vanishes identically. As the fundamental degrees of 
freedom of the theory are the partons, the presence of 
these extra bound states carrying integer electron charge 
e does not affect the physics. The important requirement 
is that Gauss' law, which enforces overall Z3 neutrality, 
ensures that the net charge of the whole sample is an 
integer multiple of the electron charge. 

For Higgs fields that transform trivially under the cen- 
ter Zat^ of SU{Nc), as is the case here, it is known 
that the confined and Higgs phases are distinct In the 
Higgs phase, the system is completely gapped even in the 
charge neutral sector because the unbroken gauge group 
Z3 is discrete. However, in the limit of infinite Higgs stiff- 
ness K — >■ 00 the system behaves like a Z^r^ gauge theory 
which has an additional gapless Coulomb phase separat- 
ing the (gapped) deconfined and confined phases for large 
enough A^c-^^'^'' It it possible that this phase persists for 
finite but large enough kM' In this case, there would be 
two distinct fTI phases with the same value of 9, one gap- 
less and the other fully gapped. We also note that the 
presence of a = tt term in the statistical Z^r^ gauge the- 
ory gives rise to the presence of oblique confined phases 
with dyon condensation in addition to the usual confined 
phase, but does not remove the gapped deconfined 
and gapless Coulomb phases corresponding to the fTI. 
Finally, one can also construct Abelian Higgs models in 
addition to the non-Abelian Higgs model discussed here, 
with [/(l)^'^"^ broken to a discrete subgroup. 



III. FRACTIONAL TOPOLOGICAL 
MAGNETOELECTRIC EFFECT 

A. Chiral anomaly 

The calculation of the effective axion angle 9 via 
the chiral or Adler-Bell-Jackiw (ABJ) anomal y^^i^^ 
was mentioned in our original work*^** following ear- 
lier discussions)^ and has been spelled out in more de- 
tail elsewhere »iii^2i2^i22r^ This argument can be used 
to obtain topological insulators in any even spacetime 
dimensioUfiiiii^iS^ but for here let us specialize to the 
case of the 3D T-invariant topological insulator. The 
goal is to calculate the contribution to the effective 9 an- 
gle of a U{1) gauge group that arises from integrating 
out a fermion of charge q. As (at least in model A) our 
partons are charged under more than one Abelian group, 
we want to calculate all terms of the form 
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where Ai is the (here, Euchdean) spacetime manifold 
and the label a runs over all the Abelian groups in the 
problem, i.e. the statistical gauge groups as well as the 
Maxwell gauge group. For example, in model A we have 
a,b G {em,A,B}. The non-Abelian case works simi- 
larly as will be discussed below. Let us denote by qf 
the charge of the ith parton under the ath gauge group 
in units of the corresponding gauge coupling (e for the 
electromagnetic sector, g for the statistical sector). In 
the Dirac kinetic term for the parton, one can write a 
complex mass term with mass M as the complex bilinear 
operator '0(ReM + i'j^lm AI)ip. The T operator takes 
M into its complex conjugate M* , so the system is only 
T-invariant if M is real. In other even spacetime dimen- 
sions it is a different discrete symmetry that takes the 
role of enforcing a real mass term.^^ Once M is real we 
see that there is a Z2 choice of mass terms: M can be 
real and positive or real and negative. The two can not 
be smoothly deformed into each other in a T-invariant 
fashion without crossing M = 0, that is, without closing 
the gap. 

As all that matters in terms of physics of interfaces is 
the difference in 9, we can choose the 0-term to be zero in 
the case that M is real and positive, which one identifies 
as the topologically trivial case. At the classical level, 
we can always perform a chiral rotation ijj — ^ e''°'''^ip to 
absorb the phase of M. This chiral rotation is a symme- 
try of the massless theory. One can think of the mass as 
a spurion, i.e. as arising from the expectation value of 
a non-dynamical background field, to restore the sym- 
metry in the massive case by letting M transform as 
M — e^'"M. If the phase of M is originally 9q, with 
6q — TT for the topologically nontrivial T-invariant in- 
sulator, M can be made real and positive by a chiral 
rotation with angle a = —9o/2. However, this chiral ro- 
tation is not a symmetry of the quantum effective action, 
as the path integral measure is not invariant. Performing 
such a rotation generates a 0-term of the form given in 
Eq. (12]), with a coefficient that is determined by a triangle 
diagram with one axial current and two U{1) currents, 

^a6-5]g^g.''^^0, (4) 

i 

i.e. it is determined entirely by the gauge group represen- 
tations to which the integrated fermions belong. For the 
non-Abelian case, the mixed 0-terms (mixed between two 
gauge groups) vanish identically due to the tracelessness 
of the representation matrices. The diagonal 0-terms are 
given by a similar formula with the charges replaced by 
the trace over the generators of the group in the repre- 
sentation of the partons. A single Dirac fermion in the 
fundamental representation contributes 9 = Oq. This cal- 
culation is robust against inclusion of interactions^! as 
recently discussed in Ref. [sl. 

From Eq. (|4]) it follows immediately that Nc partons 
of electric charge e/Nc generate a 0-term for the Maxwell 
field with 9 = 9o/Nc, that is 9 = n/Nc if the partons re- 
alize a T-invariant topological insulator. No reference to 



the statistical gauge group is necessary. In the Abelian 
case it is advantageous to ensure that no mixed 6'-terms 
involving the Maxwell field and a statistical gauge field 
are generated. Such mixed terms give rise to extra con- 
tributions to the effective 0-term for the Maxwell field 
(denoted by "em") once the statistical gauge fields are 
integrated out. The latter is a delicate thing to do in the 
deconfined phase, but these mixed terms would alter the 
topological properties of the Maxwell field. Vanishing of 
the mixed terms requires that for any a ^ em, 

Y.qtqr = 0- (5) 

i 

If all partons have the same electromagnetic charge q"" — 
1/Nc, this reads 9f ~ which is simply the require- 
ment that the electron is gauge-invariant under the sta- 
tistical gauge fields. This, by construction, is automat- 
ically satisfied in our Abelian models presented above. 
Alternatively, Eq. ^ follows from the orthogonality of 
the generators titatcm = 0, a — A, B (Sec. Ill A II) . On 
the other hand, the general non-Abelian models intro- 
duced in Sec. Ill A 31 have, for every gauge group factor 
labeled by a, a total of partons with statistical charge 
q° — and a Maxwell electric charge that is more or 

less unconstrained up to the overall condition Eq. © . In 
this case mixed terms will be generated. The final value 
quoted for 9 in Ref. 54 is only obtained after integrating 
out the statistical gauge fields in this case. 

As an explicit example, take the Abelian Nc = 3 par- 
ton model, model A (Sec. Ill A II) . Assume that all par- 
tons have a topologically nontrivial mass, i.e. ^0 — 
Using the chiral symmetry to rotate the phase of all 
three mass terms to a real and positive mass, we gen- 
erate via the ABJ anomaly an electromagnetic 0-term 
with = Ctt, where 

t^-E('^n'-^- (6) 

i 

No off-diagonal Ten A Fa or Tom A Fb terms are gen- 
erated, because the corresponding anomalies iT^lt 
and ^ - Qi'^qf vanish. In the statistical sector, nonzero 
0-terms will be generated, 

dAA - Y.(qf)'7T - Gtt, 9bb - E(gf )'7r - 6^, 

i 

Oab = Oba = qfqfTT = -3^. (7) 

These extra terms play an important role in the charge 
quantization considerations below. We note that the 
above calculations only depend on the fermion content 
of the theory, and therefore can not distinguish between 
the deconfined (model A or B) or Higgs (model C) real- 
izations of this particular fermion content. 

B. Topological quantization of the axion angle 

In a topological band insulator, the low-energy descrip- 
tion of the electromagnetic response at energy scales be- 
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low the gap is in terms of a topological field theory with 
a. Sg ^ 9E ■ B termj^ From the point of view of the low- 
energy theory, the partition function is T-invariant only 
if e*"^" (or e"'^* in Euclidean signature) is T-invariant 
which, together with the Dirac quantization of fluxes ap- 
pearing in Sg, constraints 9 to take discrete values. In 
the fTI, after integrating out the massive partons, the 
low-energy electromagnetic response still takes the form 
f?E-B, but with possibly fractional values for 9. The frac- 
tional values of 9 allowed in this case follow from mod- 
ified flux quantization conditions due to the presence of 
the statistical gauge fields. 



1. Quantization in the Abehan models 

Consider the 0-term for the Maxwell field, 
^^§^e^''^fFfj,„F\p. The quantization condition on 



^cff follows directly from the Dirac quantization con- 
dition for magnetic charges, and the argument applies 
either in Euclidean^iiS^ or Minkowski spacetime'*!'^° 
(with certain conditions on the electromagnetic fields 
in the latter case). We consider the Euclidean case 
for simplicity. With periodic boundary conditions, 
spacetime is topologically equivalent to the 4-torus T^. 
The minimal value of the spacetime integral Sg of the 

6'-term 32"^ ^^ivXpF^j^uFxp is obtained when the smallest 
allowed magnetic monopole is inserted inside both 2-tori 
and T^^ where = T^^ x T^^, and T^^ is the 
2-torus generated by directions /i and v. In Euclidean 
spacetime, all directions can be taken to be spacelike 
and thus all fluxes F^i, are magnetic. We obtain 



Sg = 



dxidx2 F12 / dx^dx^ F-^ 



34, 



(8) 



where the factor of 8 comes from the permutations of the 
e-tensorj^ If the fundamental charge is e, the smallest 
allowed magnetic monopole has magnetic flux 



27r 

F=. — = Bo, 
j'2 e 



(9) 



and Sg = i9cfi, i.e. 9cs is periodic with period 2tt. In the 
parton model however, the fundamental charge is now 
e/Nc, and the above argument would yield Jrp2 F = 
and Sg — iN^9cg, which means that 9^3 would have 
period j^. Also, the requirement that the minimal al- 
lowed magnetic Maxwell flux be NcBq and no longer Bq 
seems to be in stark contrast to the real world, where Bq 
fluxes have certainly been realized. Both of those puz- 
zles get resolved by taking into account that in models 
with fractional charges interacting with statistical gauge 
fields. Maxwell magnetic charges can be accompanied by 
"color" magnetic charges, i.e. magnetic charges of the 
statistical gauge fields . We now explain how this 



increases the periodicity of ^eff from to In order 
to do that, we must first review how the Dirac quanti- 
zation condition is modified in the presence of multiple 
Abelian gauge fields. 

For multiple Abelian gauge fields, the Dirac quantiza- 
tion condition qeqm G 27rZ with qe, qm the electric and 
magnetic fluxes, respectively (in the first example, we 
had qe = e and qm = ^), through a closed 2-manifold 
such as T^, is replaced by the more generic condition 
J2a ^elm ^ 2'kZ, whcrc again the superscript a labels 
the various Abelian gauge group factors, including the 
Maxwell gauge group. For any given gauge group a, q^qm 
does not have to be an integer multiple of 2tt. Therefore 
a minimal Bq flux is allowed, even though the product 
q^q^ of the Maxwell magnetic charge producing this flux 
and the electric charge of a parton would be 2tt/Nc, as 
long as there are also color magnetic fluxes present. This 
is possible because q^, for the parton is nonzero for the 
statistical gauge groups. Outside the topological insula- 
tor, the statistical fields are confined and their magnetic 
fiux has no physical consequence. Indeed, confinement of 
the color electric fluxes corresponds to condensation of 
the magnetic fluxes,^^ so the latter fluctuate wildly and 
there is no energy cost associated with them. Therefore, 
we recover the fact that in a topologically trivial insula- 
tor a Bq fiux is possible. Inside a fTI, the statistical U{1) 
gauge field is deconfined and the color magnetic flux does 
have an effect. 

To determine the periodicity of 0, we need to consider 
the contribution of all gauge fields to the action Sg. We 
consider the general topological term Eq. ([3|) with A4 — 
T'*, together with the result of the anomaly calculation 
Eq. (|4]). Recall that qf denotes the charge of the ith 
parton under the ath gauge group in units of the gauge 
coupling, so that q^ — eqf for the parton. Since F"' = 
q^, one finds 



a.b.i 
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The Dirac quantization condition G 27rZ en- 

sures that this is times a sum of integers squared, so 
6*0 has the standard 2tt periodicity. The Maxwell 9 an- 
gle is 9o/Nc, according to Eq. (jH), and so has periodicity 
2t: /Nc as announced earlier. 

We now demonstrate this explicitly on the example of 
model A. The Dirac quantization condition, i.e. the con- 
dition that the parton wave function should be single- 
valued, allows a Bq fiux for the Maxwell field together 
with a color magnetic flux Bq/S for the (say) U{1)b mag- 
netic field. To see this, we note that the phase by 
which the wave function of the i parton changes when 
taken around a loop enclosing this fiux is 



a.,^2n[eqf^ + eqrBQ 



(11) 
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which yields ai = 0, a2 = 27r, as = and the wave 
function is indeed single- valued. A calculation of the 9ab 
angles in model A, similar to Eq. ([7]) but with a general 
6*0, yields 6*0111 — Oa/3 and 0bb ~ 66*0 , and we obtain with 
a,b € {em, A, B} and in Lorentzian signature 



„i(6eoX^ + %xl") 



(12) 

as announced earlier, has periodicity 27r, and hence 
the Maxwell 9 angle has periodicity 2tt/3. 



2. Quantization in the non-Abelian models 

In the non-Abelian case, the discussion of charge quan- 
tization closely follows the Abelian case discussed above. 
The proper quantization condition for magnetic fluxes in 
that case is • nif, g 2TTZSab, where is a "electric flux 
vector" , nib is a "magnetic flux vector" , and the indices 
a, b run over the generators of the gauge group 
Mathematically, is a vector in the weight lattice of 
the Lie algebra of the gauge group, and the quantiza- 
tion condition defines m;, as a vector in the dual weight 
lattice, i.e. the "reciprocal" weight lattice/^ 

Let us consider the example of a deconflned non- 
Abelian gauge theory model based on TVc = 3 partons 
and an SU{3) gauge group (model B). Including electro- 
magnetism, the total gauge group is U{3) = SU{3) x 
U{l)em/^3 which has the following Cartan generators. 



H. = 
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71 ^-71 
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\ 




(13) 



where A3, Ag are Gell-Mann matrices, and we have explic- 
itly written the SU{3) and U (l)om gauge couplings g and 
e, respectively. Hi and H2 are the two Cartan generators 
of SU{3),^^ normalized to ir{HaHb) = g'^Sab, a, 6 = 1, 2, 
and H3 is the generator of ?7(l)cm with all three quarks 
having the same electric charge e/3. The weight lattice 
is generated by the fundamental weights which are 



ei 



62 



63 



9 9 e 
^/2' V6' 3 

9_ _9_ £ 

^/2' V6'3 
2.9 e 
71' 3 



0, 



(14) 



We recall that is the eigenvalue of the Cartan genera- 
tor Hi associated with the ath common eigenvector |ea) 
of all three Cartan generators, i.e. Hi\ea) = e^lea).— 
The first two entries of Bq correspond to non-Abelian 
SU (3) "color" charges, and the last entry corresponds to 
the usual t/(l)om electric charge. The dual weight lattice 
is generated by its own fundamental weights, which are 
defined as the "reciprocal lattice basis vectors" nib, 

ea ■ nib = 2TTSab, a,b ^ 1,2,3. (15) 

The linear system Eq. (fT4|) . ([T5]l is easily solved to yield^ 

/I 11' 



mi = 27r , 

\V2g VQg e 

/111 

m2 = 27r - 



V2g' V6g' e 



(16) 



The allowed magnetic monopoles, i.e. the allowed mag- 
netic flux configurations, are given by linear combina- 
tions of the fundamental dual weights mt with integer 
coefficients, 

m = ?iimi + 7i2m2 + n3m3, ni,n2,n^ G Z. (17) 

We are now in a position to discuss the periodicity of 
dcs- The "colorless" magnetic monopole configuration 
discussed earlier, which led to a periodicity of = 
corresponds to the dual weight vector 

/ 67r\ 
m = mi -f m2 -f m3 = I 0, 0, — I , 

i.e. Til = 71,2 = = 1 in Eq. (jl7p . However, we now see 
that this is not the "smallest" magnetic monopole: we 
can choose a smaller monopole for which some of the n^, 
« = 1, 2, 3 are zero. In particular, the smallest monopoles 
have only one Ui equal to 1. But as seen in Eq. 
these monopoles will necessarily be "colored", i.e. they 
will carry some amount of non-Abelian magnetic charge. 
Let us now evaluate 5*61 for a colored magnetic monopole, 
say m = mi, and see how it affects the periodicity of 
Ocs- Since the monopole is colored, we cannot simply 
discard the SU (3) 9-teim in the effective action [Eq. (6) 
in Ref.m. Denoting by F and / the C/(l)em and SU{3) 
field strengths, respectively, we have 



Sff — 
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i.e. has periodicity 27r, which means that OeS = 0/3 
has periodicity as claimed in Ref. [sj. It is readily 
checked that the other "minimal" monopoles, m = m.2 
and m = m3, give the same quantization condition. 

3. Quantization in the Higgs models 

The analysis of charge quantization in the Higgs mod- 
els is very similar to the Abelian case. The basic state- 
ment is as before. A Maxwell Bq flux is possible only 
if it is accompanied by color magnetic flux of the un- 
broken gauge group H, otherwise only NcBq is possible. 
While color magnetic flux before was quantized in inte- 
ger multiples of a basic unit, it only now takes a finite 
number of discrete values. For example, in our model C, 
magnetic flux for iJ = Z3 can only take three different 
values: three units of color magnetic flux are equivalent 
to no color magnetic flux. 

IV. GROUND STATE DEGENERACY 
A. Ground state degeneracy on 

Recently it has been shown^ that a fractional 6 an- 
gle in a T-invariant, gapped system necessarily implies 
multiple ground states on . However, a deconfined 
gauge theory may or may not be gapped, and a unique 
ground state on is possible if the system is gapless. 
In general, deconfined gauge theories with continuous 
gauge groups will be gapless and deconfined gauge theo- 
ries with discrete gauge groups will be gapped. In par- 
ticular, we will show that the Abelian models in the de- 
confined phase contain gapless but electromagnetically 
neutral gauge bosons, and have a unique ground state 
on T^. On the other hand, the non- Abelian models do 
exhibit a nontrivial ground-state degeneracy on , even 
though they avoid the theorem of Ref. [H due to the ex- 
tra massless degrees of freedom. In the Higgs models, 
the theory obeys all the assumptions of the theorem and 
there has to be multiple ground states on . We will 
see that this is indeed the case. Therefore, various re- 
alizations of a fTI with the same value of the Maxwell 
6 angle can be distinguished by their topological ground 
state degeneracy on . This is similar to the fact that 
the Hall conductance, although a topological invariant, 
is not sufficient to fully characterize the topological order 
in a FQH system. 

Before we discuss our three different models in detail, 
we must first define precisely what we mean by ground 
state degeneracy. As discussed in detail in the previous 
sections, in at least two of our three models (the models 
in the deconfined phase including models A and B), the 
theory contains gapless degrees of freedom. In order to 
arrive at a meaningful definition of topological ground 
state degeneracy in a gapless system, we need to study 
the theory for a finite-size system. Denoting by L the 



linear system size, the massless degrees of freedom will 
develop a finite-size gap of order 1/L. In order to fo- 
cus on the ground state manifold, we want to study the 
theory at energies below that finite-size gap. For a finite- 
size system, the ground state degeneracy will be generally 
lifted by nonperturbative effects^ corresponding to the 
tunneling of fractionally charged partons around noncon- 
tractible loops in , which leads to an energy splitting of 
order ~ e~™8''p^/L where mgap is the dynamically gener- 
ated parton mass gap already present at infinite volume. 
For finite L, these states are truly degenerate only in the 
strict TOgap — >■ 00 limit (i.e. for infinitely massive par- 
tons), but even at finite Wgap we can identify the ground 
states by a finite-size scaling analysis of the many-body 
spectrum. 



1. Abelian models 

We first investigate the question of ground state degen- 
eracy on in the simplest Abelian three-parton model, 
model A. When studying the model on a T'^, the ground 
state is unique. As the statistical U{1) fields are free 
in the infrared, this is simply a question of quantizing 
Maxwell electrodynamics. As is well known, C/(l)^'=^^ 
gauge theory on in the deconfined phase is sim- 
ply equivalent to a set of 2{Nc — 1) decoupled three- 
dimensional (in the x, y, z directions on T^) harmonic 
oscillators, i.e. a harmonic oscillator for each of the U{1) 
gauge fields with two polarizations. For a finite size sys- 
tem, the ground state is unique and corresponds to each 
of these harmonic oscillators being in their ground state. 
The gap to the first excited state is of order ^ 1/L where 
L is the linear system size, but vanishes in the L — ^ 00 
limit. According to Ref. |59l a gapped system can not 
have a fractional 9 angle if it has a unique ground state 
on T^. As mentioned, the way our simple Abelian model 
avoids this constraint is because it is not gapped (for 
infinite size). To see how the extra free photons can cir- 
cumvent the no-go theorem of Ref. [H, we observe that 
the argument of Ref. [5^ is essentially a Lorentzian ver- 
sion of the usual Dirac quantization argument. Consider 
a gapped scenario where a 2-torus of the spatial is 
pierced by a magnetic flux Bz- In the presence of a unit 
flux Bz = Bo = the ground-state to ground-state 
(G2G) amplitude when inserting the same minimal flux 
through the noncontractible loop in the z direction picks 
up a phase given by 

The time-reversed process picks up a phase e~*^, there- 
fore for a T-invariant theory we require e'^ = e~*^ and 9 
has to be an integer multiple of tt. 

This argument relies on the fact that the minimal mag- 
netic flux is Bq. In a theory of a Maxwell gauge field 
alone with partons of charge e/3 a magnetic fiux Bq is 
not consistent with single- valuedness of the parton wave 
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function; rather, the minimal flux allowed is SBq. In this 
case, the above G2G amplitude picks up a phase factor 
of e*^®, hence apparently any multiple of 7r/9 would be 
an allowed T-invariant value for 9. Clearly, this situation 
can not correspond to nature as we know it since a Bq 
flux can exist. At this stage, the discussion is completely 
parallel to that of Dirac quantization in Sec. IIII B II As 
we have seen in Eq. (IT^ . a basic Bq flux is allowed as long 
as it is accompanied by a magnetic flux of the statistical 
gauge fields. The Dirac quantization condition implies 
that the allowed combinations of fluxes are exactly the 
ones for which e^^ evaluates to e*"'^ for integer n. The 
G2G amplitude is T-invariant despite the fact that the 
Maxwell 9 angle takes the fractional value tt/3. In fact, 
Dirac quantization will ensure that the same happens in 
all the Abelian models proposed in Ref. 54, The no-go 
theorem is avoided even without ground state degeneracy 
on due to the presence of extra massless gauge fields. 

2. Non-Abelian models 

A deconfined SU (Nc) gauge theory on has de- 
generate ground states*2i There are two complementary 
ways of establishing this result, one by keeping the funda- 
mental partons in the spectrum, and the other by sending 
the parton mass to infinity, mgap — oo, and consider- 
ing the pure SU{Nc) gauge theory. For concreteness we 
will discuss the case Nc = 3. In the presence of funda- 
mental quarks, one can construct a topological symme- 
try algebra consisting of operators Ua, a — 1,2,3 which 
insert a 2tt fiux of the [/(l)em gauge field through the 
ath noncontractible loop of T^, and operators Ta which 
move a parton around the ath noncontractible loop. In 
the T-invariant case, we have [Ta,Tf,] = [Ua,Ub\ — 0, 
a,6 = 1,2,3, huiM 

TaUb = e-^^'/^^''UbTa, (19) 

which simply means that partons can pick up a non- 
trivial C/(l)em Aharonov-Bohm phase because they are 
fractionally charged. The operator Ta moves a parton 
around a closed loop and is clearly a symmetry of the 
partition function. The operator Ua is also a symme- 
try of the partition function because the t/(l)om phase 
g27rz/3 acquired by the partons is also an element of 
Zg = {1, e^'''/^, e'^'''/^}, the center of SU{3), and can thus 
be gauged away. The existence of an algebra of opera- 
tors which commute with the Hamiltonian but not among 
themselves implies the degeneracy of the energy eigen- 
states. The symmetry algebra Eq. (jlOp is somewhat rem- 
iniscent of the topological symmetry algebra in the FQH 
states.^*' A representation of the algebra Eq. ([T9|) can be 
constructed by first constructing a representation of the 
Abelian subalgebra generated by the Ta, i.e. by taking 
a set of states which diagonalizes the Ta simultaneously, 
Talv) — e'^^'^lv) with ij = (?7i, ?72, ?73). Pick a particular 
noncontractible loop a. By applying Ua to jrj) repeat- 
edly and using the commutation relations (I19|) . one can 
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FIG. 1: SU{Nc) gauge transformation on a noncontractible 
loop Cx in the spatial manifold [Eq. H21[)]. periodic up to an 
element 6^'''*''^'= of the center Zjv, of SU{Nc), and under 
which the global Wilson loop W{Cx) transforms nontrivially 
as W{Cx) e-^'"''^^''W{Cx). 



show that {|J7) = Ualrj) ,Ua\ri) ,Ua\ri)} is a 3-dimensional 
representation of the subalgebra generated by Ti,T2,T3 
and the Ua corresponding to this specific loop. For gen- 
eral Nc, we would obtain the A^c-dimensional represen- 
tation {|r,) = U^^\ii),Ua\rj),U^\ri), . . . , U^^-'\v)}. The 
analysis can then be repeated by starting from this en- 
larged set of states and applying the Ua corresponding 
to the remaining noncontractible loops. At the end of 
this process one finds that the dimension of the repre- 
sentation of the full symmetry algebra (|19|). which is the 

same as the ground state degeneracy, is Nc^^'^^^ where 
bi{M) = dimffi(M,M), the first Betti number^^^ of the 
spatial manifold M, corresponds to the number of non- 
contractible loops in M. It is the dimension of the first 
homology group Hi{M,M.) of M with real coefficients. 
For T^ we have bi{M) = 3. 

Alternatively, this result can be obtained by studying 
the pure SU{Nc) gauge theory which is the low-energy 
effective description at energies much less than the par- 
ton mass gap mgap. In this language, the deconfined 
phase on T^ corresponds to the condensation of spatial 
Wilson loops, i.e. a spatial version of the condensation 
of the Polyakov loop2^ (temporal Wilson loop). If all 
fields are in the adjoint (such as is the case once the 
partons have been integrated out), SU{Nc) gauge theory 
formulated on a spacetime 4-manifold Ai with nontrivial 
first homology group Hi{M) ^ develops a global 
symmetry^^'^^ originating from the center of SU{Nc). 
This global symmetry is generated by large gauge trans- 
formations, where the gauge parameter is not periodic 
on T^ but only periodic up to an element of the cen- 
ter. An order parameter for spontaneous breaking of 
this center symmetry is a spatial Wilson loop along a 
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noncontractible loop Ca, 



W{Ca) = trPexp ig 



(20) 



where g is the Yang-MiUs gauge couphng, a is the statis- 
tical SU {Nc) gauge potential, P indicates path-ordering 
along the loop, and the trace is in the fundamental repre- 
sentation. The fact that W{Ca) transforms nontrivially 
under the center of SU{Nc) is easily seen by regulariz- 
ing the theory on a lattice. Denote lattice sites by a 
triplet of integers n = (rix, Uy^Uz), with Ua = 1, . • . , Na, 
a = x^y^z, and periodic boundary conditions on the 



link variables [/„ 



Ury+Naa^i-, = t,x,y,z, in all 



spatial directions a. Consider the following family of 
Zat^ C SU{Nc) local gauge transformations pa- 
rameterized by an integer nix. 



1, 

^2TTik/Na 



(21) 



with k = 1, . . . ,Nc. Although the gauge transformation 

Eq. (PT|) is not periodic fin™"'-' ^ ^^[i+lv ^^^^^ ^ valid 
gauge transformation because the usual gauge-invariant 
operators tr[/n are invariant under any local gauge 
transformation, including multivalued ones (recall that 
the plaquette variable Ua = Un,f_,Un+ii,i.Ul_^_~^^Ul„ 

transforms in the adjoint as Utj ilnUtjfll^ for an ar- 
bitrary local gauge transformation fin)- Another way 
to say this is that although i^n"^' is not periodic as 
a SU{Nc) gauge transformation, it is periodic as a 
SU{Nc)/1'N^ gauge transformation. Since the gauge 
fields transform trivially under Z^v^ because they are in 
the adjoint, such a gauge transformation preserves the 
periodic boundary conditions on the gauge fields. 

We can now show that Eq. (pUj) transforms nontrivially 

under fin™^-*. The global Wilson loop Eq. pO| has a 
natural lattice regularization, 



W{Ca) 



tr n Ur. 



Denoting = {uy^riz), we see that W{Cx) transforms 
under fl^"'' as 

W{Cx) = tr U(i^nj_),xUi2,nj_),x ' ' ' U(N^,n^),x 

= irU(^m^^n^),xU(„i^ + l,nj_),x ' ' ' t^(m^+A'^-l,nj_) 



tr(l ■ Uf^rn^,ii^),x ■ 1)(1 ■ t^(mx + l,ni),x ' 

x(l.L^(„^^.^_i,„,).e-^-'=/^=) 



!)■ 



= e-2-«'=/^eVK(0, (23) 

using the periodicity of the trace to shift the base point of 
the loop from (l,n^) to (r7ia:,n^) and the fact that the 

link variables C/n,/j transform as Un,fj. — > ^nUn.fj.^n^fj,- 
When the theory is quantized, expectation values of 



operators such as W{Ca) are computed by averaging 
over all SU (Nc) gauge transformations that are periodic 
fin = ^n+N^a.- Note that the transformation law Eq. (|23|) 
is independent of the local data rux specifying where the 
discontinuity of the gauge function Q^iT''^ occurs. There- 
fore, although such local data can be lost by perform- 
ing the SU{Nc) gauge averaging, the global data (i.e. 
the parameter k) is not averaged out and W{Ca) can 
develop a nonzero expectation value. This is why the 
transformation law Eq. is in fact a global symme- 
try. As a result, spontaneous breaking of this symmetry 
{W{Ca)) ^ does not violate Elitzur's theoremi^ 

The Wilson loop Eq. (|20p can be interpreted as 
the semiclassical process of creating a heavy parton- 
antiparton pair at a; = and annihilating them again at 
X = L/2 along the noncontractible loop Ca, or in other 
words the worldline of a virtual parton taken around the 
loop once. Therefore, we recover essentially the same 
physics as when we explicitly considered the fundamen- 
tal partons in the previous approach of establishing the 
ground state degeneracy. The Polyakov loop is the finite- 
temperature version of Eq. (|20p . where the spatial man- 
ifold need not have a nontrivial first homology but the 
relevant noncontractible loop in that case is the periodic 
imaginary time direction r G [0,/3], with /3 = 1/T the 
inverse temperature. In this case, the logarithm of the 
Polyakov loop is the negative of the free energy asso- 
ciated with a single parton. In a confined phase, this 
free energy is infinite as it takes an infinite amount of 
energy to separate a single color charge from its color 
charge conjugate, and the Polyakov loop has zero expec- 
tation value. Center symmetry is unbroken. On the other 
hand, in a deconfined phase the Polyakov loop typically 
has a nonzero expectation value and center symmetry is 
spontaneously broken. In our case we are considering 
a zero-temperature theory and the spatial Wilson loop 
Eq. (I^H)) characterizes G2G amplitudes at zero tempera- 
ture. As each W{Ca) carries charge under Z^v^, there are 
Nc different values that W{Ca) can take for each noncon- 
tractible loop Ca- Therefore, there will be A^^ degenerate 

(22) ground states on and Nc^''^^^ ground states for a gen- 
eral spatial manifold M without boundary. 

We remark that this ground state degeneracy was not 
required by the theorem of Ref. [s^ Whether this Nc^ 
degeneracy actually occurs depends on details of the sys- 
tem beyond the number of partons and the gauge group. 
In particular, if the extra gapless neutral matter added 
in order to drive the non-Abelian statistical gauge field 
into a deconfined phase is in the fundamental represen- 
tation of the gauge group, there is no global center sym- 
metry and we expect that the ground state would be 
unique. Another example where spontaneous breaking of 
the global Zat^ symmetry (and hence presumably the cor- 
responding ground state degeneracy) does not necessarily 
occur is A/" = 4 SYM theory. As long as one imposes 
periodic boundary conditions for the adjoint fermions 
along all noncontractible loops, the potential for the cor- 
responding Wilson loop is flat due to supersymmetry. 
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The ground state expectation value of the Wilson loop is 
one more modulus in the theory. One can tune its expec- 
tation value arbitrarily, hence we can obtain a deconfined 
phase with Z^v^ either broken or unbroken. 



3. Higgs models 

Because the Higgs models of Sec. Ill Bl are completely 
gapped, they must have multiple ground states on 
to be consistent with the theorem of Ref. [H. As we 
necessarily have unbroken discrete gauge groups, this is 
indeed ensured. Discrete gauge groups give rise to degen- 
erate ground states that differ by the value of the Wilson 
loop of the discrete gauge field around noncontractible 
loops, i.e. discrete global fluxes. For example, in our 
model C with its Z3 unbroken gauge group, there are 
three different discrete fluxes per noncontractible loop of 
the spatial manifold T^, which label the various ground 
states for a total ground state degeneracy of 3"^ — 27. 
A generic SU{Nc) Higgs model on a boundaryless spa- 
tial manifold M with all fields in the adjoint, and where 
SU{Nc) is spontaneously broken to its center Zjv^ will 
have a ground state degeneracy of A^c^*-^^-*. This is the 
three-dimensional version of the N'^ ground state degen- 
eracy on in the deconfined phase of Z^r gauge theory 
in 2 -|- 1 dimensions A similar mechanism is also re- 
sponsible for the ground state degeneracy in the Z2 spin 
liquid model put forward in Ref. |59| as an explicit real- 
ization of a fully gapped fTI. 



B. Ground state degeneracy on 3-manifolds with 
boundaries 

On a 3-manifold with boundaries, such as the carte- 
sian product M = S X / of a Riemann surface S and an 
interval / discussed in our previous workj^ we need to 
also consider the effect of the Chern-Simons (CS) terms 
induced on the boundary. To be specific, let us focus on 
the M = S X / case with boundary dM = S U E. The 
nontrivial 6-term of the statistical gauge fields in the bulk 
induces a CS term on the boundary due to the axion do- 
main wall^ between 6* 7^ inside the fTI and 6 — out- 
side the fTI. For the case of 6* = tt, the corresponding CS 
term has level 1/2. The fact that this term is not gauge- 
invariant as a purely (2 -|- l)-dimensional theory does not 
matter, because the gauge noninvariance of the boundary 
theory is compensated by the bulk. This is made clear 
by writing the O-tevm as a manifestly gauge-invariant 
trt^^^Pf^^fxp term^ii22 ^^^^.^ f jg the statistical field 
strength. The statistical gauge field also has a kinetic 
(Yang-Mills) term ~ — tr/^,^/'"^. (In this section we 

work with rescaled fields — )• The details of the 

ground state degeneracy will depend on the low-energy 
dynamics of the gauge fields. We have already shown 
that the gauge fields need to be in a deconfined phase to 



have a fTI. In the following, we will further assume that 
the deconfined phase occurs at weak coupling g <^ 1, so 
that we can consider the gauge fields as being essentially 
free. The situation is less clear in the case of a deconfined 
but strongly coupled theory such as Af — 4 SYM the- 
ory. Just as in the Abelian case, the Yang-Mills Hamil- 
tonian is minimized by setting the statistical electric and 
magnetic fields to zero, which corresponds to a zero sta- 
tistical field strength = 0. Therefore, the ground 
state eigenspace is obtained by quantizing the space of 
flat connections a^, i.e. pure gauges = Ud^U~^ with 
t/ : M X K — >• SU{Nc) a smooth gauge function and R de- 
notes time. Denoting by wcs = tr (a A da -I- i|a A a A a) 
the CS form and using the fact^^ that tr / A / = dwcs 
with / = da + iaAa the statistical field strength, we have 

+ ! ua,) . (24) 

with M^_^^L EE e{z = ±L/2 + e) - e{z = ±L/2 - e), 

for e a positive infinitesimal, and denotes the CS 
form evaluated at = a^(z = ±L/2). Since z = ±L/2 
correspond to the boundary between a fTI and the vac- 
uum which is a trivial insulator, we generically have 
A0^^^ i = {2k± + l)7r where k+,k- € Z. It appears that 
we have two distinct CS terms, one for each boundary 
gauge field a^. However, those two CS terms are not in- 
dependent because the two (2-t- l)-dimensional boundary 
gauge fields originate from the same (3 -I- l)-dimensional 
gauge function U which is smooth in the bulk. As a re- 
sult, and a~ are equal up to a gauge transformation of 
even winding number which leaves the partition function 
invariant (see Eq. (|A7[) of Appendix A). The resulting 
theory is that of a single CS gauge field on S x R with 
integer level k = + + 1. This CS theory gives a 
additional contribution to the ground state degeneracy. 

In the Abelian case, we obtain a sum of U{l)k CS terms 
on the boundary (see Appendix B). However, the ground 
state degeneracy is not simply the product of the ground 
state degeneracy for each U{l)k CS term, because large 
gauge transformations which mix several U{1) factors 
give additional constraints on the ground state Hilbert 
space.''" For the Abehan [/(l)^ theory (Sec. Ill A ip cor- 
responding to Nc = 3, the ground state degeneracy is 
^{k + l)(fc + 2) for S = T2,60 which reproduces the fa- 
miliar threefold degeneracy of the = 1/3 FQH state 
on the torus for k = 1. For the Zat^ Higgs models, the 
Maxwell term of the Z^^ gauge theory itself defines a 
topological field theory and the ground state degeneracy 
is N^^ where g is the genus of E.— 

In the non-Abelian case, we obtain a non-Abelian 
CS term on The ground state degeneracy for a 

SU{Nc)k CS theory on E is equal to the number of con- 
formal blocks of the level-A: SU{Nc) Wess-Zumino-Witten 
(WZW) conformal field theoryji23, This number has been 
determined for any gauge group G, level k and genus g of 
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^ ^104,105 -j^y^-^ |.jjg answer is particularly easy for the spe- 
cial case of the torus, £1 = 1. In this case the groundstate 
degeneracy at level k is given by ^ / {Nc~l)lYl^J:^^ (k+j) . 
For the special case of fc = 1 we are mostly interested in 
this simply evaluates to Nc- For Nc = 3, the ground state 
degeneracy is ^{k + l)(k + 2) in complete agreement with 
the answer found in Ref. [60. This is the same as the 
Abehan ?7(1)^ theory. 

This discussion can be generalized to more generic 
manifolds M which are not of the product form S x /. 
Due to similar arguments as above, the ground state 
degeneracy should only depend on the "zero modes" 
with vanishing spatial gauge curvature. For each topo- 
logically nontrivial loop Ci on the surface dM, we de- 
fine a holonomyiS^ or Wilson loop operator T{Ci) = 

P exp aj . It is important to note that Ci's are con- 

sidered equivalent if they can be deformed to each other 
through some deformations in M (rather than being re- 
stricted to 9M). For example, in the case of x / the 
four cycles on the two boundary tori are identified pair- 
wise, leading to two independent nontrivial loops. After 
taking this condition into account, the operators T{Ci) 
on distinct topologically nontrivial loops have commuta- 
tion relations determined by the boundary CS theory in 
the same way as the x / case discussed above. 



V. GAPLESS SURFACE STATES 

So far we have considered effective gauge theories for 
systems with periodic boundary conditions in all spa- 
tial directions, i.e. the 3-torus T^, or for systems with 
boundaries but where the boundary is gapped since it is 
described by a CS term. In these two cases, the fermionic 
partons are gapped everywhere including on the bound- 
ary, which allows us to integrate them out. The CS terms 
break T on the boundary and are absent if T is preserved 
everywhere. In the latter case, we expect that the fTI 
should support gapless surface states since each color of 
partons condenses (at the mean-field level) into a topo- 
logical band insulator state, which does support gapless 
surface states. The question therefore arises: what is the 
nature of the gapless surface state of a fTI? 

Since we do not at present have a microscopic model of 
fTI, i.e. a model of interacting electrons, it is at present 
difficult to answer this question. Since as we have seen, 
several different effective gauge theories can give rise to 
the same quantized fractional 9 angle, we expect a variety 
of gapless surface states with properties highly dependent 
on the details of the microscopic model. 

From the point of view of the effective gauge theo- 
ries discussed in this work, the gapless surface states 
consist of a helical liquid of partons interacting with a 
three-dimensional gauge field. An effective theory for 
the (2 + l)-dimensional surface could be obtained by 
integrating out the bulk gauge fiuctuations. A similar 
calculation was performed recently^ for the surface he- 



lical spinon liquid in spin-charge separated topological 
Mott insulators)^ using a perturbative approach. In 
this case, the bulk consists of a deconfined U{1) gauge 
field. Therefore, the results of Ref. should apply 
qualitatively for the models of Sec. Ill A 11 i.e. the de- 
confined C/(l)^'=^^ models. Indeed, since the C/(l)^'=~^ 
gauge theories are deconfined at weak coupling 5 <C 1 
{g is the parton-gauge boson coupling), we expect that 
perturbation theory in g should be reliable. Further- 
more, Ref. |4^ shows that due to the three-dimensional 
nature of the gauge fluctuations, perturbation theory is 
better controlled than in the two-dimension case. How- 
ever, since the microscopic degrees of freedom in a fTI are 
electrons which correspond to gauge-invariant "baryon" 
operators in the parton gauge theory, one should only 
calculate correlation functions of gauge-invariant opera- 
tors. Ref. |4^ finds that perturbation theory at 1-loop 
gives only a logarithmic modification of the tree-level re- 
sult for the 2kp surface spin-spin correlation function, 
i.e. (S'+(r)5'_(0)) - l/ir^lnkpr). We there- 

fore expect that the 2kp surface current-current correla- 
tion function of the Abelian deconfined fTI, i.e. its sur- 
face electromagnetic response, should also only exhibit 
logarithmic modification compared to the noninteracting 
helical Fermi liquid. In the fully gapped Higgs models, 
the parton-gauge boson coupling should be irrelevant be- 
cause the bulk gauge fluctuations are massive. Therefore, 
we expect that the electromagnetic response of the gap- 
less surface state should be the same as that of the non- 
interacting helical Fermi liquid, up to corrections that 
are irrelevant at low energies. Finally, for the deconflned 
non- Abelian models, since most known examples of these 
(such as M — A SYM theory) occur at nonzero coupling 
g = g^:, we expect that the gapless surface state will be 
a strongly correlated version of the helical Fermi liquid 
and it is difficult to guess what its properties will be. We 
conjecture that the resulting state is a "helical non-Fermi 
liquid" , and holographic realizations of fTl2^ may be a 
useful tool to compute its properties. 



VI. SUMMARY 

In this work, we considered a variety of gauge theories 
in 3 + 1 dimensions and discussed the conditions they 
must fulfill to be consistent low-energy descriptions of a 
fTI. A fTI was defined phenomenologically in previous 
work^^i^ as a T-invariant state of interacting electrons 
which exhibits a quantized fractional axion angle 9 in 
its low-energy electromagnetic response. We considered 
Abelian U{1) models and non-Abelian models. In both 
cases, the confined phase is not an option for a fTI be- 
cause there would be no fractionally charged excitations 
in the spectrum, and the existence of fractionally charged 
states is necessary for a fractional 9 angle to be consistent 
with T. This leaves us with two options: a deconfined 
(or Coulomb) phase and a Higgs phase. The deconfined 
phase of Abelian U{1) models is a theory of noninteract- 
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ing, gapless gauge bosons. The gaplessness of the gauge 
bosons does not affect the quantization of 9 because they 
are electricaUy neutral. We showed the fractional quan- 
tization of 6 explicitly using the Adler-Bell-Jackiw chiral 
anomaly, which did not require the assumption that the 
gauge bosons should be gapped. Achieving a deconfined 
phase in non-Abelian models is more difficult, but non- 
Abelian models with sufficient electrically neutral gapless 
matter, such as Af = A SYM theory, are known to realize 
deconfined phases. These are however usually strongly 
coupled phases. However, the chiral anomaly still holds 
in the case of non-Abelian gauge groups, and we could 
again show explicitly the fractional quantization of 9. 
Higgs models in which the Abelian U{1) or non-Abelian 
SU{N) groups were broken down to a discrete Zjv group 
were shown to lead to fTI as well, with a fully gapped 
spectrum in this case. 

We investigated the ground state degeneracy of these 
effective gauge theories on spatial 3-manifolds of non- 
trivial topology. On the three-torus T^, the deconfined 
Abelian U{1) models have a unique ground state. In- 
deed, this corresponds simply to quantizing several in- 
dependent flavors of Maxwell electrodynamics in a box 
with periodic boundary conditions in all three directions. 
The deconfined non-Abelian models can have a nontrivial 
ground state degeneracy on due to the fact that the 
first homology group Hi{T^,Z) = Z x Z x Z is nontriv- 
ial, corresponding to the existence of three inequivalent 
noncontractible loops in T^. Whether there are multiple 
ground states or not in the deconfined phase of a SU{N) 
non-Abelian model depends on whether the center Zjv 
of SU{N) is spontaneously broken in the ground state 
or not. This is a question of dynamics which depends 
on the details of the model. Higgs models with residual 
Zjv gauge group can be viewed as a subset of the previ- 
ous case. Indeed, SU{N) gauge theories with all fields 
in the adjoint representation develop the center Z^r as a 
global symmetry which can be spontaneously broken in 
the ground state. Our Higgs models consist of adding 
adjoint Higgs fields to pure SU{N) gauge theory such 
that SU{N) is spontaneously broken to its center. Wil- 
son loops around noncontractible loops will still trans- 
form nontrivially under Zat, and their acquiring nonzero 
expectation values means spontaneous breaking of this 
global Zjv symmetry and multiple degenerate ground 
states. We restricted our consideration of 3-manifolds M 
with boundary dM to the case M = E x / with S a Rie- 
mann surface (say, in the x, y directions) and / an interval 
(say, in the z direction). In this CS term was in- 

duced on the boundary dAI = E U E, and in the ground 
state the CS gauge fields on both copies of E were iden- 
tified. The resulting CS theory had integer level and its 
contribution to the total ground state degeneracy (bulk 
and surface) could be computed using standard methods. 

Finally, we briefly commented on what one would ex- 
pect for the electromagnetic response properties of the 
gapless surface states based on the general characteris- 
tics of the effective gauge theories discussed here. We 



expect that deconflned Abelian U{1) models should give 
at most a logarithmically modifled version of the non- 
interacting helical Fermi liquid, while the fully gapped 
Higgs models should only give corrections that are irrel- 
evant at low energies. The deconfined, strongly coupled 
non-Abelian models should give rise to the most interest- 
ing case. We conjecture that the gapless surface states of 
non-Abelian fTI are "helical non- Fermi liquid" states and 
suggest that holographic methods^S, should be a promis- 
ing way to study their properties. 
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Appendix A: Flat connections and winding numbers 

In the ground state, the allowed gauge field configura- 
tions in the (3-1- l)-dimensional spacetime MxM. are pure 
gauges = Udi^JJ~^, i.e. flat connections. Denote by 
= U±dfj,U^^ the boundary values of the gauge field, 
i.e. a± = a^,{z = ±L/2) and U± = U{z = ±L/2). We 
first observe that the (2-|-l)-dimensional boundary gauge 
fields and are related by a gauge transformation, 

a+ = Loa'^^io^^ +u)d^u)^^ where lo — U+UZ^- Therefore, 
we have 

— f ujia — — [ ujro. + 2nm, (Al) 

where 

= "2^ / t^M^"')' 

is the winding number of ut. In Eq. (|24|) . the CS level 
A0^^^i/27r, which appears in front of the CS terms in 

Eq. (jAl|) . is half-integer. As a result, e''^ in the path 
integral acquires a factor of e*'^"*!'^!. If gauge transfor- 
mations uj with odd winding number m[uj] were allowed, 
the path integral would be highly oscillatory and would 
cancel. In fact, we can show that m[uj] is necessarily an 
even integer. Consider the integral of the second Chern 
form, 

/ tr/A/= / / wcs, (A3) 

JMxR JSxR JSxR 
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since dM = S U S. We now substitute the pure gauge 
a = UdU^^ into Eq. (|A3p . Since this is a pure gauge 
/ = 0, the left-hand side of Eq. (jASP vanishes. Since 



U±dUj.^, we also have f± = da± + ia± A a± = 
and the CS forms are given by 

2 2 i 

a± A da± + i-a% = a± A (/± - m|) + i-a% = ^3"±' 



hence we obtain 

i 

3 ■iT.x 







tT{U+dU-^f-^ I tT{U-dUI^f 



(A4) 



Furthermore, one can show that the winding numbers for 
a product of gauge transformations add, 



m[ujiu}2] = m[uji\ + m[uj2], 



(A5) 



as one would expect based on physical intuition. Using 
Eq. (|A2p and the cyclic property of the trace, we have 



m[uJiUJ2\ = m[wi] + m[a;2] 

+ 3/ tl [{LU2dLL>2^)'^ A duj^^LOi 

+u}2duj2^ A (dwf . (A6) 

Using the identity dr] = —ridrj^^rj for any invertible 
matrix-valued function -q, we can show that 

d {LU2dLjJ2^ A duJi^UJl) = — {ijJ2dijJ2^Y ^ duj^^LUi 
— UJ2dljJ2^ A [dljj^^ LOi)^ , 

hence the last term of Eq. (|A6p is a total derivative and 
Eq. (|A5p is obtained. Also, using the same identity one 
can easily show that m[ri^^] = ~m[rj\. For cj ~ U+UZ^ 
we therefore obtain 

m[(jj\ — m[J7+] + to[[/~"'"] ~ to[[/_|_] — m[J7_], 

and from Eq. (jA4p we have m[C/+] + m[?7_] — 0, hence 

m[uj] = 2m[U+\ = -2to[C/_] e 2Z. (A7) 



Appendix B: Ground state degeneracy for an 
Abelian 3D fTI on a 3-manifold with boundary 

We consider the Abehan C/(l) x t/(l) fTI on the spatial 
3-manifold M = x I with / = [~L/2,L/2] in the z 
direction. Because the bulk 3D statistical gauge fields are 
free, we can integrate them out explicitly to obtain an 
effective action for the gauge fields on the 2D boundary 
dM = U T^. It is a Maxwell-CS theory with two 
coupled gauge fields a+ and a~ corresponding to the two 
copies of . In the long- wavelength limit q ^ 1/L the 
two gauge fields become identified = = and 
the level of the CS term for is the sum of that for the 
two surfaces, i.e. it is integer. 



The spatial manifold is M = x / with / = 
[— i/2,L/2]. We can always choose the generators of 
U{1) X U{1) to satisfy trtah = Sab, a,b = 1,2. The 
action in imaginary time is 



1 



ie{z) 



ra ra 

2Jti,vJti.v 327r2 "^^ 



with Z;^^ = - d^a"^, a = 1,2 the U{1) x U{1) "sta- 
tistical" field strength and 

d,9 = (2^ + - rjL/2), 

n=±i 

with fc,, e Z. To derive an effective 2D action on 
dM ~ T^UT^, we introduce a Lagrange multiplier which 
constrains the gauge field to live on dM. Then we inte- 
grate out the bulk gauge field a°. In other words, we 
introduce a resolution of unity 



1 = J V~a+V~a-5[al{i,L/2) - h^-^i)] 

x5[al{i,^L/2)-~a-^%i)], 

in the partition function, where aJJ'°, 77 = ± are two 
auxiliary gauge fields which are defined only on the 2D 
surface x — (xq — t,xi,X2), with a^'° living on the 2- 
torus at z = L/2 and a"'" living on the 2-torus at z = 
—L/2. We represent the functional delta function as 

X exp [i J d^ij;'^ii)[a;{i,rjL/2) - d^'^ii)]^ , 

where is a Lagrange multiplier which implements 
the constraint that d^''^{x) = aj^(i, 77L/2). The idea is to 
integrate out first a^, and then the Lagrange multiplier 
to get an effective action 821) in terms of the 2D gauge 
fields d'^ alone. In other words, the partition function is 



Z = / T^%T^d^ Vj^Vj^ Va^ exp {^-Ssd [a^] 

+ . / d^xj;^'^{x)[a;{i,nL/2)-di^'^{i)] 



We want to see whether these gauge fields will be iden- 
tified or not, i.e. whether we obtain two CS theories or 
only one. 

First of all, since the 6'-term is a total derivative it 
contributes only to the boundary piece. Therefore if we 
write = + S'jq we immediately know that 

'-'2D — * 2-^ ^ / " xe^y\a^ o^a^ , 
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i.e. we obtain two decoupled CS terms of half-odd- integer 
level kjj + ^. To compute S^^^ we need to integrate 
out the 3D bulk gauge fluctuations, which we can do 
exactly in the Abelian case because the gauge bosons are 
noninteracting. We can gauge-fix the Maxwell term in 
the usual way'^'^^ by adding a 2p{(^A''^/i)^ term to the 
Lagrangian. In the Feynman gauge ^ = 1 the 3D gauge 
boson propagator is 

D-l{k,h)^l^^, (Bl) 

with k — (fco — w,fci,fc2). Because the Lagrange multi- 
pliers jji'"'{k) = J (Px e~'^'^"^'' j^'°'{x) are independent of 
fcz, we have to sum over all kz to obtain an effective 2D 
propagator. The effective 2D Maxwell propagator for the 
2D gauge fields aJJ'° is therefore 

/ ~ ,1 p-ikAv~v')L/2 

with = v^cj^ + k2. Because the fields aJJ^'^ must satisfy 



some sort of boundary condition at z = ±L/2 (either 
Dirichlet or Neumann), k^ is a discrete variable, fc^ = 
nir/L, n e Z. Performing the discrete sum over kz, we 
obtain 

r^ah ( cosh\q\L 1 \ 
^''^^> \q\smh\q\L\ 1 coshlqlL^' 

where we denote q = k for simplicity. The inverse prop- 
agator is 

^ ( cosh|g|L -1 \ 
\^ -1 cosh|g|L j ■ 

To obtain the final form of 82%, we need to integrate 
out the Lagrange multipliers j^'" which simply amounts 
to inverting the 2x2 matrix propagator (jBip . Since 
now all quantities are 2D, we can drop all the tildes for 
simplicity, and obtain 



^2DK,«^J-2„2r / ^27r)3sinhlalL/'"'rA. )- 



2g-^Lz J {2^fsiii\i\q\Lz \^ ^ J -<i\ -I co8h\q\Lz J \ a''^ J ^ 

lY^hll fd^,,^^,al'^d.ar, (B2) 



with q= (w, q) and \q\ — \/ uj"^ + q^. We denoted Lz = L for clarity. 

Consider fixing the scaling dimension of the gauge fields by the CS term. The latter is therefore marginal and 
contains one power of q. The effective Maxwell term in Eq. (jB2[) contains higher powers of \q\. Let us expand it to 
quadratic order in 

^eff r^-f ,-1 ^ f ( a-- 1 ( ' + + ^'^'^^ ^ ( ^ 



We now consider the long-wavelength, low-energy limit 
\q\ <Si 1/ Lz- In this limit, the quadratic Maxwell terms 
q^L^ are irrelevant and the leading term is 




2g^L 



This term contains no derivatives of aj^'° but simply im- 
plements a constraint. The equations of motion read 



rcoff 
rccff 

Salj.' 



dL 
dL 

da.- 
flu 



2«'+-«r)' 



-2(a^'+-a'l--), 



which imply that a"'"*" = a^' and the gauge fields on the 
two 2-tori are identified. Therefore, in the limit \q\ <C 
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l/Lj the CS term in Eq. becomes with a° = a°'+ = a°~ and = + = k+ + k- + l 

is the effective CS ievel, which is integer. Therefore, the 



■'2D 



E'"') + 5 /" 73 jj a Q v,a ground state degeneracy of the U{1) x U{1) Abehan fTI 

47r " J 2;ep^AV '^'^^a' on x J is the same as that of the level-A; J7(1) x [/(I) CS 

theory ([BSl on T^, which is known^ tobei(A:+l)(fc+2). 
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